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shown that embedding diagrams are closely related to the notion of the radius of 
gyration which is useful for analyzing fluid rotating in strong gravitational fields. 
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1. Introduction 

The optical reference geometry related to stationary spacetimes enables to introduce 
the concept of inertial forces in the framework of general relativity in a natural way 
[1, 2]. Of course, in accord with the spirit of general relativity, alternative approaches 
to the concept of inertial forces are possible (see, e.g., [3, 4]); however, here we shall 
follow the approach of Abramowicz and his coworkers [5], providing a description of 
relativistic dynamics in accord with Newtonian intuition. 

The optical geometry results from an appropriate conformal (3 + 1) splitting, 
reflecting some hidden properties of the spacetimes under consideration through its 
geodesic structure. The geodesies of the optical geometry related to static spacetimes 
coincide with trajectories of light, thus being 'optically straight' [6, 7]. Moreover, 
the geodesies are 'dynamically straight,' because test particles moving along them are 
kept by a velocity-independent force [8]; they are also 'inertially straight,' because 
gyroscopes carried along them do not precess along the direction of the motion [9] . 

Some properties of the optical geometry can be appropriately demonstrated 
by embedding diagrams of its representative sections [1, 10, 11]. Because we are 
familiar to the Euclidean space, usually 2-dimensional sections of the optical space 
are embedded into the 3-dimcnsional Euclidean space. (Of course, embeddings into 
other conveniently chosen spaces can also provide interesting information, however, 
we shall focus our attention on the most straightforward Euclidean case.) In the 
Kerr-Newman backgrounds, the most representative section is the equatorial plane, 
which is their symmetry plane. This plane is also of great astrophysical importance, 
especially in connection to the theory of accretion disks [12]. 

In the spherically symmetric spacetimes (Schwarzschild [6], Reissner-Nordstrom 
[10], and Schwarzschild-de Sitter [13]), an interesting coincidence appears: the turning 
points of the central-plane embedding diagrams of the optical space and the photon 
circular orbits are located at the same radii, where, moreover, the centrifugal force, 
related to the optical space, vanishes and reverses sign. 

However, in the rotating black-hole and naked-singularity backgrounds, the 
centrifugal force does not vanish at the radii of photon circular orbits in the equatorial 
plane [14]. Of course, the same statement is true if these rotating backgrounds carry 
a nonzero electric charge. It is, therefore, interesting to study how the inertial forces, 
defined in the framework of the optical geometry, and the photon circular orbits are 
related to equatorial-plane embedding diagrams of the optical geometry of rotating, 
charged backgrounds. Such relations were discussed in the case of non-charged, Kerr 
backgrounds in [15]. In this paper, we shall generalise the results to the more complex 
case of the Kerr-Newman backgrounds, in which even stable photon circular orbits 
can exist beside the unstable ones, contrary to the case of Kerr backgrounds [16]. 

In Section 2, the optical reference geometry and the related inertial forces are 
defined relative to the family of locally non-rotating observers in the Kerr-Newman 
spacetimes. In Section 3, stationary equatorial circular motion in the Kerr-Newman 
spacetimes is discussed, using the concept of the gravitational and inertial forces 
expressed in terms of a 'Newtonian' velocity related to the optical geometry. It is 
shown that asymptotically the relativistic expressions of the gravitational, Coriolis 
and centrifugal forces reduce to the well known Newtonian formulae. In central parts, 
they enable an illumination of unusual properties of the Kerr-Newman spacetimes in 
terms of intuitively clear concepts. The properties of the centrifugal force in the 
equatorial plane are closely related to the properties of the embedding diagrams 
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of the optical-geometry equatorial plane. In Section 4, the embedding formula 
is introduced, the limits of Euclidean embeddability of the optical geometry are 
established, and the turning points of the embedding diagrams are determined; it is 
also shown that they occur just where the centrifugal force reverses sign. The locations 
of photon circular orbits in the equatorial plane are given, and it is established 
which of them arc contained in the embeddable regions of the optical space. The 
Kerr-Newman spacetimes are classified according to the criterion of embeddability 
of regions containing photon circular orbits. Finally, typical embedding diagrams 
are constructed, and the Kerr-Newman spacetimes are classified according to the 
properties of the embedding diagrams (namely, the numbers of embeddable regions 
and turning points of the diagrams). In Section 5, some concluding remarks are 
presented. 



2. Optical geometry and inertial forces 

The notions of the optical reference geometry and related inertial forces are convenient 
for spacetimes with symmetries, especially for stationary (static) and axisymmetric 
(spherically symmetric) ones. However, they can be introduced for a general spacetime 
lacking any symmetry [2]. 



2.1. General case 

Assuming a hypersurface globally orthogonal to a timelike unit vector field n K and a 
scalar function <!> satisfying the conditions 

n[ K V\n^ = 0, n K n K = -1, h x = n K V K n x = V A $, (1) 

the 4- velocity u K of a test particle of rest mass m can be uniquely decomposed as 

u K = 1 {n K + vr K ). (2) 

Here r K is a unit vector orthogonal to n K , v is the speed and 7 = (1 — u 2 )~ 1//2 . 

Introducing, according to Abramowicz, Nurowski and Wex [2], a projected 3-space 
orthogonal to n K with the positive definite metric giving so called ordinary projected 
geometry 

h K \= g K \ + n K n x , (3) 

and the optical geometry h K \ by conformal rescaling 

h K \ = c~ 2 *(g KA + n K n x ), (4) 

the projection of the 4-acceleration aj^ = u^^ux can be uniquely decomposed 
into terms proportional to zeroth, first and second powers of v, respectively, and the 
velocity change 

i = (cV),/. (5) 

Thus, we arrive to covariant definition of inertial forces analogous to the Newtonian 
physics [2, 17] 

maj: = G K (v°) + C^v 1 ) + Z K {v 2 ) + E K (v), (6) 
where the first term 

G K = -mV B $ = -m$ K (7) 
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corresponds to the gravitational force, the second term 

C K = -m 7 2 vn A (V A T K - V K r A ) (8) 
corresponds to the Coriolis-Lense-Thirring force, the third term 

Z K = -m( 1V ) 2 T A V A f K (9) 
corresponds to the centrifugal force, and the last term 

E K = -mvf K (10) 

corresponds to the Euler force. Here f K is the unit vector along r K in the optical 
geometry, and V K is the covariant derivative with respect to the optical geometry. 

2.2. Kerr-Newman case 

Using geometric system of units (c = G = 1), and denoting M the mass, a the 
specific angular momentum, e the electric charge, the line element of the Kerr-Newman 
spacetime, expressed in terms of standard Boyer-Lindquist coordinates, reads 

As 2 = - f 1 J At 2 i — '- sin 2 9 dtdct> + 

^|^d^ + |dr 2 + Sd^ (11) 

where 

A = r 2 -2Mr + a 2 + e 2 , (12) 

E = r 2 + a 2 cos 2 6, (13) 

A=(r 2 + a 2 ) 2 -Aa 2 sin 2 d. (14) 

For simplicity we put M = 1 in the following; equivalently, we use mass units of 
M . If a 2 + e 2 < 1, the metric (11) represents black-hole spacetimes. The loci of 
their horizons, r_ (the inner one) and r + (the outer one), determined by real roots of 
A(r; a, e) — 0, can be equivalently given by the relation 

a 2 = al(r;e) = r(2 - r) - e 2 . (15) 

The case a 2 + e 2 = 1 corresponds to extreme black holes. If a 2 + e 2 > 1, there are no 
horizons, and the metric (11) represents a naked-singularity spacetime. If e 2 < 1, both 
the black-hole and naked-singularity spacetimes contain an ergosphere, where gu < 0; 
particles and photons in bound states with covariant energy E < are possible there 
[18]. Naked-singularity spacetimes with e 2 > 1 have no ergosphere. 

The Kerr-Newman spacetimes, being stationary and axially symmetric, admit 
two commuting Killing vector fields; the vector field r\ K is (at least asymptotically) 
time-like, having open trajectories, the vector field £ K is spacelike, having closed 
trajectories. Now, the vector field n K relevant for constructions of the ordinary 
projected geometry and the optical reference geometry can be given by using these 
Killing vector fields [5] , and corresponds to the family of locally non-rotating frames 
(LNRF) or zero angular momentum observers (ZAMO) introduced by Bardeen [12]. 
Namely, 

n K = C*(77 K + Q L NR F e), ^LNRF = -eVx/V^, /, fi v 

$ = -\ in (-7?V - 20 LNRF £V - nl NRF e*i\) ■ { 1 
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The LNRF vector field n K can be used for the definition of inertial forces introduced 
above. Assuming a circular motion with angular velocity = d(f>/dt as measured by 
the stationary observers at infinity, the 4-velocity is given by 

u^i^ + OD, A= (-r?V-2fi£V-nV&)" 1/2 ; (17) 

now t k is directed along the rotational Killing vector £ K . The gravitational (7), 
Coriolis-Lense-Thirring (8), and centrifugal (9) forces can be written down as 

9t<t> ~ 9tt94>4> \ 



In . 

\ 9<t><t> J 



(18) 



C K = mA 2 [(etx)WZJ, K - (rf&ietxU = 

mA 2 (tt - ^lnrf)V500 d K (gt<j>9^ 2> j + ^lnrf d K ^/g^ , (19) 



Z K . — 



^mA 2 (n - fi L NRF) 2 30* d K 




(20) 

9ttg<t>4> I 

respectively; we denote i K = e~*n K . The Eulcr force will appear for £1 ^ const 
only, being determined by Ct — u x \7\tt (see [5]). The electromagnetic forces acting in 
charged spacetimes were defined in the framework of the optical geometry in [17, 19]. 
However, here we shall concentrate on the inertial forces acting on the motion in 
the equatorial plane (9 = 7r/2), and their relation to the embedding diagram of the 
equatorial plane of the optical geometry. 

Clearly, by definition, the gravitational force is independent of the orbiting 
particle's velocity. On the other hand, both the Coriolis-Lense-Thirring and the 
centrifugal force vanish (at any r) if Cl = £1 LNRF , i.e., if the orbiting particle is 
stationary at the LNRF located at the radius of circular orbit; this fact clearly 
illustrates that the LNRF are properly chosen for the definition of the optical geometry 
and inertial forces in accord with Newtonian intuition. Moreover, both the forces 
vanish at some radii independently of fi. In the case of centrifugal force, namely this 
property will be imprinted into the structure of the embedding diagrams. 

3. Stationary equatorial circular motion in the Kerr Newman spacetimes 

For the stationary circular motion, O = 0. It is convenient to express the inertial 
forces in terms of 'Newtonian' velocity 

v = "fv. (21) 

Velocity v takes values from — oo to oo, while v from —1 to 1. In the stationary and 
axially symmetric spacetimes, there is 

v = QR, (22) 

where 

= 0-0 LNRF , R = ?e*, f=(C^ K ) 1/2 - (23) 

R is so called radius of gyration, because R 2 = £/Q, where £ = L/E is the specific 
angular momentum, L = t/ K £ K is the angular momentum, and E = —U K r} R — 7e* 
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is the energy [5, 7]. It plays a very important role in theory of rotational effects in 
strong gravitational fields. The direction of increase of the radius of gyration gives 
a preferred determination of the local outward direction relevant for the dynamical 
effects of rotation. This direction becomes misaligned with the 'global' outward 
direction in strong fields [7] . The condition R = const defines the von Zeipel cylinders 
in stationary spacetimes, which are related to the equipotential surfaces of equilibrium 
configurations of perfect fluid [20]. 

It is convenient to introduce the gravitational acceleration, and the velocity 
independent parts of the Coriolis and centrifugal accelerations in the direction e K 
by the relations [5] 

g(r)=e K V K $, C(r) = e K .RV K LNRF , Z(r) = e K R~ 1 W K R. (24) 

The acceleration necessary to keep a particle in a stationary motion with a velocity v 
along a circle r = const in the equatorial plane can then be expressed in a very simple 
form 

a(v, r) = -g(r) - v 2 Z(r) + (1 + v 2 ) 1/2 iC{r) (25) 

that enables an effective discussion of the properties of both accelerated and geodesic 
motion. (Of course, only the positive root of the last term on the r.h.s. of Eq. (25) 
has physical meaning.) 

For the stationary circular motion in the equatorial plane of the Kerr-Newman 
spacetimes all three parts of the acceleration have only radial components. We obtain 

_ r 4 (r - e 2 ) + a 2 [2r(r - 2){r - e 2 ) - e 4 ] + a 4 (r - e 2 ) 
y[T) ~ rA[(r 2 +a 2 ) 2 - a 2 A] ' ( ' 

= 2a[r(3r 2 +a 2 )-e 2 (2r + a 2 )} 
1 ' r\/A[(r 2 +a 2 ) 2 - a 2 A] ' 1 ' 

Z(r) = {rA[(r 2 + a 2 ) 2 - a 2 A}}' 1 {r\r 2 - 3r + 2e 2 )+ 

a 2 [r 2 (r 2 - 3r + 6) + e 2 r(3r - 7) + 2e 4 ] - 2a 4 (r - e 2 )} . (28) 

These definitions of the gravitational and inertial forces have really a Newtonian 
character, since the gravitational force G = —mQ(r) is velocity independent, C = 
m(l + w 2 ) 1 / 2 wC(r) depends on v, and the centrifugal force Z = —mv 2 Z(r) depends 
on v 2 . Moreover, asymptotic behaviour of these forces is consistent with 'Newtonian' 
intuition: 

G(r ->oo) 4. C(r ->oo) 4-, Z(r -» oo) ~ -. (29) 

It follows immediately from Eq. (25) that the photon circular geodesic motion 
(v 2 — > oo) is determined by the conditions 

Z(r) — C(r) = 0, corotating orbits, (30) 

Z(r) + C(r) = 0, counterrotating orbits. (31) 

For ultrarelativistic particles (v 3> I , v <§C — I ), we obtain an asymptotic relation 

a(v, r) w -Q(r) ± \C(r) - v 2 [Z(r) T C(r)]. (32) 

The upper signs correspond to the corotating motion (y > 0), the lower signs to the 
counterrotating (v < 0) motion. The ultrarelativistic particles moving on the radius of 
the corotating (counterrotating) photon circular geodesic are kept by the acceleration 

a{r) = -g(r)±\C{r), (33) 
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which is independent of velocity, with accuracy 0(v 2 ). In static spacetimes C(r) = 0, 
and at the radius of the photon circular orbit the acceleration of particles is 

a(r) = -G(r), (34) 

and is independent of v exactly. (The relation (34) is not limited to the case of 
ultrarelativistic orbits, because Z(r) = at the radius of the photon circular geodesies 
in static spacetimes.) 

Velocities of particles moving along the circular geodesies are determined by the 
relation 



C 2 -ZGt \C{C 2 - 4ZG + AG 2 ) 1 ' 2 



*± = ' 2 ' Z 2_ C 2 , (35) 

where the upper (lower) signs correspond to the corotating (counterrotating) orbits, 
v > (v < 0). 

Properties of the gravitational acceleration Q(r), and of the velocity independent 
parts of the Coriolis and centrifugal of acceleration C(r) and Z{r) determine properties 
of both the accelerated, and geodesic motion in a given spacetime. It is clear that for 
the geodesic circular motion at a given r, the condition 

C 2 (r) - AZ{r)G(r) + 4£ 2 (r) > (36) 

must be satisfied. In static spacetimes C(r) = 0, and the reality condition of the 
geodesic orbits is Q(r)Z(r) < 0, i.e., the gravitational and centrifugal forces must 
point in opposite directions. If C(r) ^ 0, the cases \Z\ > \C\, \Z\ = \C\, \Z\ < \C\ 
give qualitatively different situations. Of course, the gravitational force can play 
an important role too. Therefore, we have to determine behaviour of the functions 
Q(r;a,e), C(r;a,e), and Z(r;a, e). We give an illustration of the behaviour of these 
functions for a black-hole (Fig. la) and naked-singularity (Fig. lb) spacetimes. The 
asymptotic behaviour of these functions is given by the relations (29). Further, it is 
important to know, if these functions are positive or negative. 

We have to find, where the functions Q(r;a, e), C(r;a,e), and Z(r;a,e) change 
their sign. All of these functions diverge, and change their sign, at the boundary of 
the region of causality violations, which is determined by the relations 

a 2 = ag. v .(r;e) = e2 _^_ r2 . (37) 




Figure 1. Behaviour of the velocity-independent parts of gravitational (solid), 
Coriolis (dashed), and centrifugal (dashed-dotted) forces Q(r), C(r), and Z(r), 
respectively, (a) a 2 = 0.16, e 2 = 0.16, (b) a 2 = 0.86, e 2 = 0.16. 
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In the Kerr spacetimes (e 2 = 0), this region is restricted to r < 0. 

Surprisingly, the gravitational acceleration changes its sign at the zero points 
given by 



a 2 = al ± (r; e) = v ' v 6 , (38) 



_ e 4 - 2r(r - 2)(r - e 2 ) ± ^/D~ g 
^ [nc '- 2(r - e 2 ) 



with 



£> g = 4r(r - 2)(r - e 2 )(2re 2 - 2r 2 - e 4 ) + e 8 . (39) 

It holds at r > even for the Kerr spacetimes — we arrive at a simple formula 

a 2 = a 2 g± {r)=r(2-r± 7) . (40) 

The Coriolis acceleration changes sign at zero points determined by 

2 2/ x ^(3r 2 -2e 2 ) 
a 2 ^ Q g(r; e )^ \ 2 _ r ■ (41) 

In the Kerr spacetimes, C(r) does not change sign at r > 0. 

Finally, we find that the centrifugal acceleration changes sign at radii given by 

a 2 = at ± (r; e) = ^ - 3r + g) + 4 ^- ^ + ± ^, (42) 

where the discriminant is 

D z (r; e) = r 8 + 2r 7 - (3 + 2e 2 )r 6 + 4(2e 2 - 9)r 5 + 
3(12 + 26e 2 -e 4 )r 4 - 

2e 2 (42 + 27e 2 )r 3 + e 4 (73 + 12e 2 )r 2 - 28e 6 r + 4e 8 . (43) 

We give an example of the behaviour of all the functions a 2 ± (r;e), a 2 (r;e), and 
a 2 ± (r; e) in Fig. 2. Generally, properties of the stationary circular motion can be given 
in a straightforward manner by using these functions; we shall not discuss details here. 

The centrifugal acceleration is very important, since it is closely related to the 
radius of gyration (see Eq. (24)). Therefore, it deserves a detailed study. 




Figure 2. Behaviour of the functions a^ ± (solid), (dashed), and a^ ± (dashed- 
dotted). They are given for fixed charge parameter of the spacetime e 2 = 0.16. 
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Figure 3. The function determining reality condition D z (r;e) > for the 
radii at which the centrifugal force vanishes. Inside the shaded region, there 
is D z (r;e) < 0, and the functions a 2 ± (r;e) are not denned there. The three 
remaining solutions of the equation D z (r;e) = are not presented here, being 
physically irrelevant (with a 2 < or e 2 < 0). 

The reality condition for a 2 ± (r; e) is 

D z (r;e)>0, (44) 

which can be treated easily considering D z (r; e) as a quartic polynom in e 2 . The results 
are given in Fig. 3. Detailed behaviour of the functions a 2 ± (r;e) will be determined 
in the next section, since it is closely related to embedding diagrams of the equatorial 
plane of the optical geometry. 




4. Embedding diagrams of the optical geometry 

Metric coefficients of the optical geometry of the Kerr-Newman spacetimes are, due 
to Eq. (4), given by the relations 

I hee A ~ A 2 . 2 

h rr = — = ^ , h H = sin 9. (45) 
In the equatorial plane they reduce to 

h rr (9 = ^tt) 



(r 2 +a 2 ) 2 -a 2 A 



A 2 



h*W = ¥) - dm -■ ( 46 ) 



_ ^ _ [(r 2 + a 2 ) 2 -a 2 A] : 
r'A 

The properties of the optical geometry related to the Kerr-Newman spacetimes 
can conveniently be represented by embedding of the equatorial (symmetry) plane 
into the 3-dimensional Euclidean space with line element expressed in the cylindrical 
coordinates (p, z,<f>) in the form 

dcr 2 = dp 2 + p 2 d<P + dz 2 . (47) 
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The embedding diagram is characterised by the embedding formula z — z(p) 
determining a surface in the Euclidean space with the line element 

dp 2 + p 2 dtf (48) 



d£ 2 E) 



dp J 



isometric to the 2-dimensional equatorial plane of the optical space determined by the 
line element 

d£ 2 = h rr dr 2 + h u dcj) 2 (49) 

with h rr and given by (46). 

The azimuthal coordinates can be identified, the radial coordinates are related as 

p 2 = hfa, (50) 

and the embedding formula is governed by the relation 

dZ X = hrr(^) 2 -l. (51) 



^dp J \dp 

It is convenient to transfer the embedding formula into a parametric form z{p) 
z(r(p)) with r being the parameter. Then 




(r 2 + a 2 ) 2 - a 2 A ( dp 



A 2 Kdr < • < 52 > 

The sign in this formula is irrelevant, leading to isometric surfaces. Because 
dz dz dr 
dp dr dp ' 

the turning points of the embedding diagram, giving its throats and bellies, are 
determined by the condition dp/dr = 0, where 

= { r 4[ r(r _ 3 ) + 2e 2 ] + a 2 [2e 4 + e 2 r(3r - 7) + r 2 (r 2 - 3r + 6)] - 

dr 

2a 4 (r - e 2 )}[r 3 (r 2 - 2r + a 2 + e 2 ) 3/2 ]. (54) 

By comparing Eqs (28) and (54) we can immediately see that turning points of the 
embedding diagrams are really located at the radii where the centrifugal force vanishes 
and changes sign. Thus, we can conclude that just this property of embeddings of the 
optical geometry of the vacuum spherically symmetric spacetimes (see [1, 10, 13]) 
survives in the Kerr-Newman spacetimes. However, photon circular orbits are 
displaced from the radii corresponding to the turning points of the embedding 
diagrams. Therefore, it is interesting to find the situations where the photon circular 
orbits lie within the regions of embeddability of the optical geometry. 

4-1. Limits of embeddability 

The embeddability condition (dz/dr) 2 > (see (52)) implies the relation 
E(r; a, e) = 4r n - 3(e 2 + 3)r 10 + 12(a 2 + e 2 )r 9 - 

2(17a 2 + 5a 2 e 2 + 2e 4 )r 8 + 4(9a 2 + 3a 4 + 12a 2 e 2 )r 7 - 
(33a 4 + 66a 2 e 2 + lla 4 e 2 + 17a 2 e 4 )r 6 + 
4(9a 4 + a 6 + 13a 4 e 2 + 10a 2 e 4 )r 5 - 
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(36a 4 + 12a 6 + 78a 4 e 2 + 4a 6 e 2 + 21a 4 e 4 + 8a 2 e 6 )r 4 + 
2(12a 6 + 42a 4 e 2 + 12a 6 e 2 + 27a 4 e 4 )r 3 - 
(4a 8 + 52a 6 e 2 + 73a 4 e 4 + 12a 6 e 4 + 12a 4 e 6 )r 2 + 
4(2a 8 e 2 + 9a 6 e 4 + 7a 4 e 6 )r - 

4(a 8 e 4 + 2a 6 e 6 + a 4 e 8 ) > 0. (55) 

The function E(r;a,e) can be considered as a polynom quartic both in a 2 and e 2 . 
For e 2 = 0, the function E(r;a) is still a polynom quartic in a 2 (see [15] for details). 
However, for a 2 = 0, the function E(r; e) simplifies significantly, being only a polynom 
quadratic in e 2 

E(r; e) = r 8 [r 2 (4r - 9) + 3r(4 - r)e 2 - 4e 4 ] ; (56) 

its behaviour is discussed in [10]. For a 2 = e 2 = we arrive at the well-known 
Schwarzschild condition r > |. 

The limits of embeddability are determined by the condition E(r; a, e) = which 
can be solved as a quartic equation in a 2 . Let us denote the four solutions as a 2 fc (r; e), 
k e {1,2,3,4}. Instead of giving long explicit expressions for the four solutions 
a 2 fc (r;e) we will treat them numerically and classify qualitatively; different types of 
their behaviour will be described. Of course, we naturally restrict our attention to 
physically relevant situations when a 2 > and e 2 > 0. 

Recall that if e 2 = 0, the limits of embeddability are given by the solutions a 2 3 (r), 
a 2 4 (r) representing two branches, both reaching zero at r — (see Fig. 4). The upper 
branch (a 2 4 ) diverges for r — > oo. The lower one (a 2 3 ) has a local minimum at r = 1, 
where a 2 ( min ) = 1, and two local maxima, a 2 ^ max 2 ^ = 1.13540 at r = 0.68947, and 
a e(max i) = 1-07543 at r = 1.33172; its second zero point is located at r — 2.25, 
corresponding to the Schwarzschild case [15]. 

If e 2 > 0, the limits of embeddability still consist of two branches. We can consider 
three qualitatively different situations. 




Figure 4. Square of the specific angular momentum of the Kerr backgrounds 
corresponding to the event horizons (a^(r;e), the border between the dark gray 
dynamic area, in which the optical geometry is not defined, and the light gray 
non-embeddable area), photon circular orbits (Op h (r;e), dashed curve), turning 
points of the embedding diagram (a 2 ± (r;e), solid curve), and embeddability 
region border (a 2 (r; e), dotted curve enclosing the light gray area not embeddablc 
into Euclidean space) are drawn as functions of the radius for fixed value e 2 = 0. 
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4-1.1. Class Ea: e 2 E. (0,0.68950) The two branches are given by the solutions 
a 2 3 (r; e) and a 2 4 (r; e), again. Both the branches diverge at r = e 2 . The upper branch 
(a 2 4 ) diverges also for r — > 00, having a local minimum a 2 4( . min -j near r = e 2 . The 
lower branch (a 2 3 ) has always a local minimum at r = 1, where a 2 ( min ) = 1 — e 2 , 
corresponding to extreme black-hole states, and a local maximum a 2 ( max y > 1 — e 2 
at r > 1. Its zero point, corresponding to the Rcissncr-Nordstrom case, is determined 
by E(r\ e) = (cf. Eq. (56)). It can also have a local minimum a 2 3 ( min ) and a local 
maximum a 2 3( - max 2) at r < 1. Positions of a 2 4 ( min ) and a 2 3 ( min ) relative to the value 
a 2 = 1 — e 2 , and the relation between a 2 ,, . s and a 2 ,, „>, or a 2 Q , ,n, determine 

' c4(mm) c3(max,2)' c3(max,l)' 

different members of embeddable regions of the Kerr-Newman black-hole and naked- 
singularity spacetimes. We can give the following subclassification of the spacetimes 
listed by the number of embeddable regions according to the values of the parameter 
e 2 ; for fixed e 2 from given interval, the numbers of the embeddable regions are given 
with the parameter a 2 growing: 



Subclass 


Interval of e 2 


Black holes 


Naked singularities 


Eai 


(0,0.17230) 


1,2,3 


4,3,2 


Ea 2 


(0.17230,0.17906) 


1,2 


3,4,3,2 


Ea 3 


(0.17906,0.19749) 


1,2 


3,2,3,2 


Ea 4 


(0.19749,0.20670) 


1 


2 ^ 3^ 2 ^ 3} 2 


Ea 5 


(0.20670,0.20794) 


1 


2,1,2,3,2 


Ea 6 


(0.20794,0.28961) 


1 


2,1,2,1,2 


Ea 7 


(0.28961,0.68950) 


1 


2,1,2 



Notice that in all of the subclasses Eai-Ea7, both the branches of the limit of 
embeddability are defined at r > e 2 . The behaviour of the embeddability limits 
in the subclasses Eai, Ea3, and Ea6 is illustrated in the Figs 5, 6, and 7, respectively. 

4.1.2. Class Eh: e 2 <E (0.68950,2.205) The limits of embeddability have two 
branches. The upper one is determined by the solution a 2 4 (r;e). Again, it is well 
defined at r > e 2 only, diverges at r — e 2 and for r — > 00, having a minimum near 
r = e 2 . The lower branch is radically different from the class Ea. The solution a 2 3 (r; e) 
diverges at r — e 2 again, but it has a discontinuity, which is 'filled up' by various 
combinations of the solutions a 2 1 (r;e), a 2 2 (r;e), and a 2 4 (r;e). The combinations 
depend on the parameter e 2 . However, it is not worth to discuss them explicitly 
because they have common basic properties. They have no local extrema, but two 
lobes — an external (upper) one, and an internal (lower) one. The internal lobe can 
enter the region of r < e 2 for e 2 high enough; for e 2 > 1, the internal lobe is shifted 
to the physically irrelevant region, where a 2 < 0. If e 2 < 1, the solution a 2 3 (r;e) has 
a local minimum at r = 1, with a 2 = 1 — e 2 corresponding to an extreme black-hole, 
and a local maximum a 2 ( max u > 1 — e 2 at r > 1; for e 2 = 1 these local extrema 
coalesce at r = 1, a 2 3 = 0. Therefore, the subclassification according to the number 
of embeddable regions can be given in the following simple way (with a 2 growing): 



Subclass 


Interval of e 2 


Black holes 


Naked singularities 


Ebi 


(0.68950,1) 


1 


2,1,2 


Eb 2 


(1,2.205) 


none 


1,2 



An example of subclass Ebi can be seen in Fig. 8. 
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0.5 1 1.5 2 2.5 3 3.5 



r 

Figure 5. Square of the specific angular momentum of the Kerr-Newman 
backgrounds corresponding to the event horizons (a 2 (r; e), between the dark gray 
dynamic area, in which the optical geometry is not defined, and the light gray non- 
embeddable area), photon circular orbits (Op h (r; e), dashed curve), turning points 
of the embedding diagram (a 2 ± (r;e), solid curve), and embeddability region 
border (a 2 (r;e), dotted curve enclosing the light gray area not embeddable into 
Euclidean space) are drawn as functions of the radius for fixed value e 2 = 0.16. 
This case covers all four types of black-hole embedding diagrams and the first three 
cases of naked-singularity diagrams. They are determined by the lines a 2 = const, 
and depicted by the notation of the classification according to the properties of 
embedding diagrams (BH1-BH4, NS1-NS3). The dotted vertical line is common 
vertical asymptote of a 2 h (r;e), a 2 ± (r;e), and a 2 (r;e). 




r 



Figure 6. The functions a 2 (r; e), a ph (r; e), a 2 ± (V; e), and a 2 (r; e) drawn for fixed 
value e 2 = 0.19. The naked-singularity classes NS4-NS6 are depicted. 




r 



Figure 7. The functions a 2 (r; e), a ph (r; e), a 2 ± (r; e), and a 2 (r; e) drawn for fixed 
value e 2 = 0.25. The naked-singularity classes NS7-NS10 are depicted. 
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Figure 8. The functions a 2 (r; e), <J 2 h (r; e), a 2 ± (r; e), and a 2 (r; e) drawn for fixed 

value e 2 = 0.92. Contrary to the preceding cases, there exist lobes of a 2 (r;e), 
and a 2 ± (r;e) is defined for all r > e 2 . However, there exist an interval of the 
parameter a 2 corresponding to naked-singularity spacetimes, at which no turning 
points appear on the embedding diagrams — these are the NSn-spacetimes. 

4-1.3. Class Ec: e 2 6 (2.205, oo) The limits of embeddability a 2 (r;e) have two 
branches, again. However, now both the branches are determined by the solution 
a 2 4 (r;e). The first branch is defined at r > e 2 , diverges at r = e 2 and for r — > oo, 
having a local minimum near r — e 2 . The second branch is relevant from some r < e 2 , 
where a 2 4 = 0, and it diverges at r — e 2 , having no local extreme. The subclassification 
according to the number of embeddable regions is simply given just by one case (with 
a 2 growing): 



Subclass 


Interval of e 2 


Black holes Naked singularities 


Eci 


(2.205, oo) 


none 1 , 2 



Note that in situations which arc necessary in order to construct typical 
embedding diagrams, the curves a 2 (r;e) are given explicitly (see Figs 5-8). The 
behaviour of the embeddability limits of the class Ec is not illustrated because they 
give no qualitatively different kind of embedding diagrams. 

4-2. Turning points of the embedding diagrams 

The embeddable regions are characterised by the number of radii where the embedding 
diagrams have turning points, or, equivalently, where the centrifugal force determined 
by (28) vanishes. Therefore, the turning points of the diagrams are governed by 





NSn 
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the functions a 2 ± (r;e) determined by (42) and (43). We shall discuss the behaviour 
of these functions and introduce a corresponding classification of the Kerr-Newman 
backgrounds (relative to their parameter e 2 ) according to the number of turning points 
of the embedding diagrams of their optical geometry. 

It follows from the reality condition (44) (see also Fig. 3) that for e 2 < e 2 r = 
0.51033 the functions a 2 ± (r;e) are not defined between the radii n(e), r (e), i.e., 
inside the shaded region in Fig. 3. If e 2 = (see Fig. 4), there is a 2 ± (r = 0) = 0, 
and a 2 _(r = 3) = 0. The part of a 2 _(r, 0) starting at r = is located in the region 
between the horizons, and, therefore, is physically irrelevant; a local maximum of 
a%+(r) is located at r( maX; i) = 0.81159 with a 2 +( - max ^ = 1.36668; a 2 +( r ) an d a l-( r ) 
coincide at r = 1 with a 2 = 1, corresponding to an extreme black-hole (see [15]). 

For e 2 > 0, we can separate qualitatively different classes of the behaviour of 
a 2 ± (r; e) in the following way. 

4.2.1. Class Ta: e 2 6 (0,0.51033) According to the reality condition (44), the 
functions a 2 ± (r;e) are defined at r < n(e), and r > r Q (e). For r = r;(e), r = r Q (e), 
there is a 2 + = a 2 _. At r > r (e), the function a 2 „(r;e) has a zero point, while 
a 2 +(f ; e) diverges for r — > +00. The functions a 2 -(-(r; e) have no extrema at r > r (e). 
At r < n(e), the function a 2 _(r;e) grows from its zero point up to r = n(e), it 
is physically irrelevant up to r = 1, where a 2 _(r = l,e) = 1 — e 2 . The function 
a 2 + (r;e) diverges at r = e 2 ; it has a local minimum o 2 + ^ min ^(e) at r z ( m ; n )(e) and a 
local maximum a 2 +(max) (e) at r c(max) (e); of course r z(min) < r z(max) < r ; . Relating 

a z+(min) ( e ) to a 2 = 1 — e 2 enables us to give the classification of black-hole and naked- 
singularity backgrounds according to the number of turning points of the embedding 
diagrams. (A classification of these backgrounds relating the number of turning- 
points and the number of embeddable regions can be realized by relating a 2 4( min ) 

With fl c+(max). and a e3(max,2) with fl z+(min)-) 

The subclassification according to the number of turning points can be given in 
the following way (with a 2 growing): 



Subclass 


Interval of e 2 


Black holes 


Naked singularities 


Tai 


(0,0.18275) 


1,3 


4,2 


Ta 2 


(0.18275,0.51033) 


1 


2,4,2 



Examples of subclasses Tai and Ta 2 are given in Figs 5 and 6 or 7, respectively. 



4.2.2. Class Tb: e 2 6 (0.51033,1.125) If e 2 equals to the critical value e 2 r = 
0.51033, n(e) = r (e), and the two branches of a 2 + and a 2 _ coalesce. For e 2 > e 2 r , 
the reality condition (44) is satisfied everywhere, and a 2 _(r; e), and a 2 + (r; e) give two 
separated branches of the turning points. The function a 2 _(r; e) determining the lower 
branch is physically irrelevant at r < 1, being located between the horizons, and it 
has a local maximum a z_( max )(e) at r > 1. The function a 2 + (r;e) giving the upper 
branch diverges at r — e 2 and for r — > 00, and it can have a local minimum a 2 + ( min ) 
near r — e 2 , or two local minima a 2 + ^ min ^, a 2 + ( min 2 ) an d a local maximum a 2 + ^ max ^ 
between them. Now, the subclassification according to the number of turning points 
can be given as follows (with a 2 growing): 
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Subclass 


Interval of e 2 


Black holes 


Naked singularities 


Tbj 


(0.51033,0.52480) 


1 


2,0,2,4,2 


Tb 2 


(0.52480,1) 


1 


2,0,2 


Tb 3 


(1,1.125) 


none 


2,0,2 



An example of subclass Tb 2 is given in Fig. 8. 



4-2.3. Class Tc: e 2 £ (1.125, oo) Now, the function a 2 _(r;e) is irrelevant, being 
shifted to negative values completely. The function a 2 + (r; e) behaves in the same way 
as that of the subclass Tb 3 . The classification according to the number of turning 
point is the following: 



Subclass 


Interval of e 2 


Black holes Naked singularities 


Tci 


(1.125, oo) 


none 0, 2 



We can conclude that outside the outer black-hole horizon, there is always an 
embeddablc region, covering exterior of the black hole except a small part in the 
vicinity of the outer horizon, and containing just one turning point corresponding to a 
throat of the embedding diagram. So, the situation is the same as for Schwarzschild, 
Reissner-Nordstrom, and Kerr black holes. On the other hand, under the inner 
horizon, the embeddable region contains no turning point, or two turning points 
corresponding to a throat and a belly; we do not consider the case of coalescing 
of the two turning points in an inficx point separately. 

For the naked-singularity backgrounds, the situation is much more complex, 
including various possibilities of the number of the embeddable regions and the turning 
points of the diagrams. There are a lot of cases that cannot appear with Kerr naked 
singularities [15]. As examples, let us mention one region with four turning points, 
and one region with no turning point. Again, in situations which are necessary in 
order to determine behaviour of typical embedding diagrams, the curves a 2 ± (r; e) are 
given explicitly (see Figs 5-8). 



4-3. Embeddability of photon circular orbits 

The motion of a photon in the equatorial plane of Kerr-Newman spacetimes is 
determined by the function (see, e.g., [16]) 

R(r; a, e, E, $) = [Er 2 - a($ - aE)} 2 - A($ - aE) 2 , (57) 

where E is the covariant energy of the photon, and $ is its axial angular momentum. 
For photon circular orbits the conditions 

an 

R = 0, ^ = (58) 

must be satisfied simultaneously. The equatorial photon motion is fully governed by 
the impact parameter 

$ 

£=-. (59) 

It follows from the conditions (58) that the radii of photon circular orbits are 
determined by the equation 

r 2 - 3r + 2a 2 + 2e 2 ± 2aA 1/2 = 0; (60) 
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(61) 



the corresponding impact parameter is given by the relation 
_ a(r 2 + 3r - 2e 2 ) 
r 2 — 3r + 2e 2 

Equivalently, the radii of photon circular orbits can be determined by 

2 2 , x (r 2 -3r + 2e 2 ) 2 
a = a ph (r;e) = 4(r _ g2) • (62) 

Clearly, the circular orbits must be located at 

r > e 2 . (63) 
The zero points of a 2 h (r; e) are located at radii 

r 1 (e)-i[3-(9-8e 2 ) 1 / 2 ], r 2 (e) = ±[3 + (9 - Se 2 ) 1 / 2 ], (64) 

giving photon circular orbits of Reissner-Nordstrom spacetimes. The extrema of 
a 2 h (r;e) are located at the zero points r\(e), r 2 (e) (if e 2 < |), at r = 1 (if e 2 < 1) 
and at 

r = \e\ (65) 



3" 

where a minimum exists for e 2 < |, a maximum for | < e 2 < |, and a minimum for 
e 2 > |. At r = |e 2 , the value of a 2 h (r; e) is given by the function 

<&(«0 = i e 2 (8e 2 -9) 2 . ( 66 ) 
The function a 2 x (e) determines boundary between Kerr-Newman spacetimes 
containing different numbers of photon circular orbits. If e 2 < |, there can be 2 
or 4 circular orbits in the black-hole spacetimes and 2 orbits in naked-singularity 
spacetimes. If | < e 2 < 1, there are 2 circular orbits in black- hole and 2 or 4 orbits in 
naked-singularity spacetimes. For the naked-singularity spacetimes with 1 < e 2 < |, 
there are 2 or 4 orbits; if e 2 > |, there are or 2 orbits. 

In the case of extreme black holes (a 2 + e 2 = 1), the radii of photon circular orbits 
and corresponding values of the impact parameter are 

r = 1, £ = a + l/a, 

r = 2 - 2a, i = 4 - 3a, (67) 
r = 2 + 2a, £=-(4 + 3a), 

The counterrotating orbit at r = 2 + 2a is always located above the horizon. The 
corotating orbit at r = 2 — 2a is located above the horizon if e 2 > |, and under the 
horizon if e 2 < |. There are two orbits at r = 1 if e 2 < |, but no circular orbit at 
r = 1, if e 2 > | (see [16] for details). 

Now, we shall discuss in which cases the photon circular orbits enter the regions 
of embeddability. Radii of the photon circular orbits are determined by the function 
a 2 h (r;e) given by Eq. (62). The embeddability of these orbits must be determined 
by a numerical procedure. Note that, generally, the function a 2 _(r;e) has common 
points with a 2 h (r; e) at its zero points (if e 2 < |), and at r = 1 (if e 2 < 1). 

The classification of the Kerr-Newman spacetimes according to the embeddability 
of the photon circular orbits can be again given in a simple way. In the presented 
classification scheme related to the parameter e 2 , the first digit gives the number 
of circular orbits, the digits in parentheses determine changes of the number of 
embeddable circular orbits, with parameter a 2 growing. Note that the circular orbits 
under the inner horizon are always outside the embeddable regions: 
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Figure 9. Embcddability of photon circular orbits in the embedding diagrams 
of the optical geometry. It is determined by the relations of the curves a 2 (r;e) 
(dotted lines) giving the limits of embcddability, and a 2 h (r;e) (dashed lines) 
giving radii of the photon circular orbits. We present all of the typical situations 
representing the classes Pa-Pj, restricting our attention on the relevant parts of 
the functions a 2 (r; e) and i 2 h (r; e). Note that the circular orbits under the inner 
black-hole horizon are always outside the embcddable regions. 
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0.185 0.190 0.195 0.200 

a 2 



Figure 10. Classification of the Kerr-Newman spacetimes according to embed- 
dability of photon circular orbits. The solid lines separate regions of the parameter 
a 2 -e 2 corresponding to spacetimes containing different numbers of the circular 
orbits; the bold solid straight line correspond to extreme black-hole spacetimes. 
The dashed line separates regions of the parameter space with different numbers of 
the embeddable circular orbits. The S-shapcd solid curve is given by Eq. (66). The 
corresponding region of the parameter space are depicted by pairs of digits, the 
first one expressing number of the circular orbits, the second one (in parentheses) 
the number of the embeddable orbits, (a) the global view of the classification, 
(b,c) details of some narrow regions located in the vicinity of the extreme black- 
hole states, where the structure is rather complicated. 



Class 


Interval of e A 


Black holes 


Naked singularities 


Pa 


(0,0.23083) 


2(2), 4(2,1) 


2(1) 


Pb 


(0.23083,|) 


2(2,1),4(1) 


2(1) 


Pc 


(|,0.80115) 


2(2,1) 


4(1), 2(1) 


Pd 


(0.80115,0.80797) 


2(2) 


4(2,1), 2(1) 


Pc 


(0.80797,0.82040) 


2(2) 


4(2, 3), 2(1) 


Pf 


(0.82040,0.82488) 


2(2) 


4(4, 3), 2(1) 


Pg 


(0.82488, 1) 


2(2) 


4(4), 2(2,1) 


Ph 


(1 £ ) 


none 


4(4), 2(2,1) 


Pi 


(§,1.70233) 
(1.70233, oo) 


none 


0(0), 2(2,1) 


Pj 


none 


0(0), 2(2) 
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The embeddability of the photon circular orbits can be easily read out from the 
sequence of figures (Fig. 9) representing the classification given above. 

Details of the classification, i.e., distribution in the parameter space of the Kerr- 
-Newman backgrounds, must be determined by a numerical procedure. The results 
are presented in Fig. 10. The regions of the parameter plane a 2 -e 2 are denoted by two 
digits. The first one gives the number of photon circular orbits in the corresponding 
background, the second one (in parentheses) gives the number of embeddable orbits. 

4-4- Construction of the embedding diagrams 

The relevant properties of the embedding diagrams are determined by the functions 
a 2 (r;e) and a 2 (r;e) governing embeddable parts of the optical reference geometry 
and turning points of the diagrams. We shall present a classification of the Kerr- 
-Newman spacetimes according to the number of the embeddable regions, and the 
number of their turning points. All cases of the classification will be represented by a 
typical embedding diagram. Inspecting all of the types of the behaviour of functions 
a 2 ± (r; e) (given by classes Ta-Tc and their subclasses), we find that there are 4 cases 
of the behaviour of the embeddings for the black-hole spacetimes, and 11 cases for the 
naked-singularity spacetimes. A complete list of typical embeddings will be given by 
using behaviour of the characteristic functions a 2 (r;e) and a 2 (r;e) for some typical 
values of the parameter e 2 — see Figs 5-8, where the functions a\{r;e) and a^ h (r\e) 
are included for completeness. In order to obtain all of the 15 types of the embedding 
diagrams, it is necessary to consider at least four values of the parameter e 2 belonging 
subsequently to the subclasses Ear (e 2 = 0.16), Ea 3 (e 2 = 0.19), Ea 6 (e 2 = 0.25), 
Ebi (e 2 = 0.92) of the classificat ion according to the number of embeddable regions, 
discussed above. The behaviour of o 2 (r;e) for these four values of e 2 enables us to 
present a complete list of typical embedding diagrams. The classes will be denoted 
successively with growing parameter a 2 for each fixed e 2 . 

In the regions of the optical geometry where the embeddability condition 
E(r; a,e) > is satisfied, the embedding diagram can be constructed for a fixed 
parameter a by integrating the paramctrically expressed embedding formula z(r) (cf. 
Eq. (52)), and transferring it into the final form z(p) by an appropriate numerical 
procedure using (50). In all the typical diagrams also photon circular orbits are 
illustrated, if they enter the embeddable regions. (Of course, the presented diagrams 
do not cover all of the possibilities for the embeddability of photon circular orbits.) 

All of the black-hole diagrams are determined by Fig. 5. In order to cover all of 
the naked-singularity diagrams all Figs 5-8 have to be used. The typical diagrams 
arc represented in Figs 11 and 12 for the black- hole spacetimes, and in Figs 13 and 
14 for the naked-singularity spacetimes. We do not explicitly consider the cases 
corresponding to situations when extrcma of a 2 (r;e) coincide in an inflcx point as 
they can be discussed in a straightforward manner. 

The embeddability of photon circular orbits can be systematically discussed by 
using Figs 9 and 10, and behaviour of a 2 (r; e) and a^ h (r; e) in related situations. We 
shall drop off the discussion here, we only explicitly mention two interesting cases of 
naked-singularity spacetimes containing four photon circular orbits all of which are 
located in the embeddable parts of the optical geometry instead. (Note that such a 
situation cannot occur in Kerr backgrounds.) In the first case (a 2 = 0.10005, e 2 = 0.9), 
there are two separated parts of the diagram, both containing two of the circular orbits 
(see Fig. 15), in the second case (a 2 = 0.105, e 2 = 0.9), there is only one region of 



Optical reference geometry of the Kerr-Newman spacetimes 



21 




Figure 11. Embedding diagram of the Kerr-Newman black holes of the type 
BHi, constructed for a 2 = 0.16, e 2 = 0.16. The rings in the 3D diagram 
represent photon circular orbits. Both corotating (gray spot in 2D diagram) and 
counterrotating (black spot in 2D diagram) are displaced from the throat of the 
diagram, where the centrifugal force vanishes. This is a general property of the 
rotating backgrounds. 
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Figure 13. Embedding diagram of the NSi type, constructed for a 2 = 0.86, 
e 2 = 0.16. 



embeddability containing all of the circular orbits (see Fig. 16). In both cases, the 
embedding diagrams necessarily have a throat and a belly. 

4-5. The classification 

The following classification will be made according to the properties of the embedding 
diagrams. The Kerr-Newman spacetimes are characterised by the number of 
embeddable regions of the optical geometry, and by the number of turning points 
at these regions, given successively for the regions with descending radial coordinate 
of the geometry (which are presented in parentheses). Therefore, the classification can 
be given in the following way: 



BHs 




NSs 




NSs 




NSs 




BHi: 


1(1) 


NSi: 


4(1,1,1,1) 


NS 5 : 


3(1,1,2) 


NS 9 : 


2(2,0) 


BH 2 : 


2(1,0) 


NS 2 : 


3(2,1,1) 


NS 6 : 


2(2,2) 


NSi : 


1(4) 


BH 3 : 


2(1,2) 


NS 3 : 


2(3,1) 


NS 7 : 


2(1,1) 


NS n : 


1(0) 


BH 4 : 


3(1,1,1) 


NS 4 : 


3(1,1,0) 


NS 8 : 


1(2) 







The parameter space of the Kerr-Newman spacetimes can be separated into 
regions corresponding to the classes BH1-BH4, and NSi-NSn by a numerical code 
expressing the extreme points of the functions a^(r;e) and a^(r;e) as functions of 
e 2 . Results of the numerical code are presented in Fig. 17, which represents the 
classification completely. 

5. Concluding remarks 

In the rotating black- hole and naked-singularity backgrounds, the embedding diagrams 
of the optical reference geometry reflect immediately just one important property of 
the inertial forces related to a test-particle circular motion. Namely, the turning 
points of the embedding diagrams occur just at radii where the centrifugal force 



Optical reference geometry of the Kerr-Newman spacetimes 23 




(f) (g) (h) 0) (j) 

Figure 14. Embedding diagrams of the types (a) NS2 (a 2 = 1, e 2 = 0.16), 
(b) NS 3 (a 2 = 1.05, e 2 = 0.16), (c) NS 4 (a 2 = 0.82, e 2 = 0.19), (d) NS 5 
(a 2 = 0.847, e 2 = 0.19), (e) NS 6 (a 2 = 0.9, e 2 = 0.19), (f) NS 7 (a 2 = 0.77, 
e 2 = 0.25), (g) NS 8 (a 2 = 0.9, e 2 = 0.25), (h) NS 9 (a 2 = 0.99, e 2 = 0.25), 
(i) NS10 (a 2 = 1-23, e 2 = 0.25), and (j) NSn (a 2 = 1, e 2 = 0.9). 

vanishes and reverses sign, independently on the velocity of the motion. However, 
contrary to the spherically symmetric black-hole and naked-singularity backgrounds, 
in the rotating backgrounds radii of the photon circular orbits do not coincide with 
the radii of vanishing centrifugal force, and, moreover, in a variety of Kerr-Newman 
backgrounds, some of these orbits even are not located in the the regions of the optical 
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Figure 15. Embedding diagrams of the NS7 type, constructed for a 2 = 0.10005, 
e 2 = 0.9. The diagram is separated into two regions having one throat and one 
belly. All four photon circular orbits find themselves in the embeddable region. In 
each region, one corotating orbit (gray spots) and one counterrotating one (black 
spots) exist. 

geometry that are embeddable into the 3-dimensional Euclidean space. 

The properties of the embedding diagrams have been discussed, and a classifica- 
tion scheme of the Kerr-Newman backgrounds reflecting the number of embeddable 
regions, and the number of turning points of these regions, was presented. The 
ring singularity in all of the rotating backgrounds, and the horizons of the black- 
hole backgrounds, are always located outside the embeddable regions. Further, 
embeddability of regions containing the photon circular orbit has been established. 

The presence of a nonzero charge parameter e in the Kerr-Newman backgrounds 
enriches significantly the variety of the embedding diagrams in comparison with the 
case of pure Kerr backgrounds [15]. 

For the black-hole backgrounds, there is just one embeddable region outside the 
outer horizon, containing just one throat. However, under the inner horizon, there 
can be two, one, or no embeddable regions. They can have a throat and a belly. 
Above the outer horizon, the outermost, counterrotating photon circular orbit is 
always embeddable, while the inner, corotating one is embeddable in backgrounds 
with sufficiently small specific angular momentum, but non-embeddable in the other 
ones. On the other hand, the photon circular orbits existing under the inner horizon 
are always located in the non-embeddable regions. 

For the naked-singularity backgrounds, the variety of possible embedding dia- 
grams is much more complex and can be directly read out from Figs 13-14. Let us 
point out the most important phenomena. There can exist diagrams consisting from 
four separated regions, or, oppositely, simply connected diagrams, having two throats 
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Figure 16. Embedding diagrams of the NSs type, constructed for a 2 = 0.105, 
e 2 = 0.9. The diagram is simply connected having one throat and belly. All four 
photon circular orbits find themselves in the embeddable region. The corotating 
and counterrotating circular photon orbits are displayed as gray and black spots, 
respectively. 

and two bellies. On the other hand, there are also diagrams having no turning point. 
Moreover, in some naked-singularity backgrounds containing four photon circular 
orbits, all of the orbits are located in the embeddable regions — such situation is 
impossible in the black-hole backgrounds. 

The embedding diagrams of the optical geometry give an important tool of 
visualisation and clarifying of the dynamical behaviour of test particles moving 
along equatorial circular orbits, by imagining that the motion is constrained on 
the surface z(p) [10]. The shape of the surface z(p) is directly related to the 
centrifugal acceleration. Within the rising portions of the embedding diagram, the 
centrifugal acceleration points towards increasing values of r, and the dynamics of test 
particles has essentially Newtonian character. However, in the descending portions of 
the embedding diagrams, the centrifugal acceleration has a radically non-Newtonian 
character, as it points towards decreasing values of r. Such a kind of behaviour appears 
where the diagrams have a throat or a belly. At the turning points of the diagram, 
the centrifugal acceleration vanishes and changes its sign, i.e., dr/dp = dz/dp = 0. 

We can understand this connection of the centrifugal force and the embedding of 
the optical space in terms of the radius of gyration representing rotational properties 
of rigid bodies. In Newtonian mechanics, it is defined as the radius R of the circular 
orbit on which a point-like particle having the same mass M and angular velocity 
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Figure 17. Classification of the Kerr— Newman spacetimes according to prop- 
erties of the embedding diagrams. The solid lines separate spacetimes with 
different numbers of embcddablc regions, the dashed lines separate spacetimes 
with different numbers of turning points (throats and bellies, in other words), 
at which the centrifugal force vanishes. The parameter plane a 2 -e 2 is divided 
into regions corresponding to types BH1-BH4 and NSi-NSn. The central part 
marked by dotted rectangle is zoomed in the lower graph. 
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as a rigid body would have the same angular momentum J: 
J = MR 2 n. 

Defining a specific angular momentum I = J/M, we obtain 



(68) 



R = 




(69) 



For a point-like particle moving on r = const, there is R = r. 

In general relativity, the radius of a circle can be given by two standard ways — 
namely as the circumferential radius, and as the proper radial distance. A third way 
can be given by generalising Eq. (69) for a point particle moving along r = const; we 
use I = L/E, where L is the angular momentum of the particle, and E is its energy. [In 
stationary spacetimes, the angular velocity have to be related to the family of locally 
non- rotating observers, cf. Eq. (23).] In Newtonian theory, these three definitions 
of radius issue identical results, however, in general relativity they are distinct. The 
radius of gyration R is convenient for discussing the dynamical effects of rotation; the 
direction of increase of R defines local outward direction of these effects. The surfaces 
R = const, called von Zeipel cylinders, were proved to be a very useful concept in the 
theory of rotating fluids in stationary, axially symmetric spacetimes. In Newtonian 
theory, they are ordinary straight cylinders, but their shape is deformed by general- 
relativistic effects, and their topology may be non-cylindrical. There is a critical family 
of self-crossing von Zeipel surfaces [20] . 

It is crucial that in the Kerr-Newman spacetimes 



i.e., the embedding diagrams of the equatorial plane of the optical geometry are 
expressed in terms of the radius of gyration. Since in the equatorial plane the 
centrifugal acceleration has the radial component 



the relation of the embedding diagrams, the radius of gyration, and the centrifugal 
force is clear. Note that the turning points of the embedding diagrams determine both 
the radii where the centrifugal force changes sign, and the radii of cusps where the 
critical von Zeipel surfaces are self-crossing. Therefore, the embedding diagrams also 
reflect properties of fluid rotating in Kerr-Newman spacetimes. 

Notice that the black-hole backgrounds have an unified character above the event 
horizon — there exists a throat of the embedding diagram indicating change of sign 
of the centrifugal force nearby the event horizon. On the other hand, the naked- 
singularity spacetimes give a wide variety of the behaviour of the embeddings and 
centrifugal forces, ranging from the simple 'Newtonian' backgrounds with no change 
of sign of centrifugal force, to very complicated backgrounds, where the sign is changed 
four times. 

We established also the embeddability of circular photon geodesies. There is a 
lot of case when these orbits are located in regions that are non-embeddable. 

Notice that because the photon circular geodesies are given by the condition 
Z(r) — C{r) = for corotating orbits, and Z(r) + C(r) = for counterrotating orbits, 
the corotating (counterrotating) orbits are located on the descending (rising) portion of 
the embedding diagrams (assuming C{r) < 0), if they enter the embeddable regions — 
see Figs 15 and 16. 



[ft^(0 = V2)] 1/2 = P = R, 



(70) 



Z(r) = R~ 1 d r R, 



(71) 
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Finally, it should be noted that in the case of the spherically symmetric 
spacetimes (Schwarzschild and Reissner-Nordstrom [10, 18, 21], and Schwarzschild- 
-de Sitter [13, 22]), it can be shown that the limits of embeddability of the optical 
geometry coincide with the existence limit of static configurations of uniform density. 
Unfortunately, we haven't found any evident connection between the embeddability 
limits of the diagrams and any other phenomenon in the Kerr-Newman spacetimes. 
Search for such phenomenon remains a big challenge for future investigation. 
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